Quantumlike Diffusion over Discrete Sets 
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Abstract 

In the present paper, a discrete differential calculus is introduced and used to describe dynam- 
ical systems over arbitrary graphs. The discretization of space and time allows the derivation of 
Heisenberg-like uncertainty inequalities and of a Schrodinger-like equation of motion, without need 
of any quantization procedure. 
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Ever since the seminal work of Regge on gravity, 
discretization has been often used in classical physics to 
explore complex continuous theories in a noperturbative 
fashion. The idea there is to reduce the infinite number of 
degrees of freedom of (typically Riemannian) manifolds 
by dealing instead with piecelinear spaces described by 
a finite number of parameters. For example, in standard 
Regge calculus a manifold is approximated by a simpli- 
cial lattice with fixed coordination number, while in the 
dynamical triangulated random surface method i3| the 
very coordination numbers are treated as dynamical de- 
grees of freedom. The notion of discrete ambient space 
(or space-time) has thus become one of the standard ways 
of making complex geometries and topologies accessible 
to (essentially combinatorial) discrete calculations, in a 
setting naturally renormalized. 

On the other hand, when quantum features are consid- 
ered, discreteness plays an unexpected role. It was shown 
that merely importing conventional quantum models on 
ambient spaces that are represented by inhomogeneous 
graphs 13, the graph inhomogeneity plays the role 
of curvature, giving rise to effective interactions among 
free particles. The question raises of course of whether 
or not a quantum model of this sort of a system living 
on a graph, represented e.g. in the second quantization 
formalism, actually describes the real (possibly interact- 
ing) physical system we intend to represent; whether or 
not is has a classical counterpart from which it can be 
obtained by a process of quantization. Canonical quanti- 
zation, however, requires an underlying structure that is 
not manifestly evident when the system lives on a graph. 

In this paper we intend to explore some of the struc- 
tural features that underly the problem. 

Let us consider a dynamical system whose ambient 
space, discrete, is provided by a graph X, either finite 
or infinite. The graph will be assumed to be completely 
defined by its vertex and edge sets, V{X) and E{X), re- 
spectively. The edge connecting two given vertices v and 
v' will be denoted as {vv'). The geometry of X arises 



completely from the adjacency relations encoded in the 
adjacency matrix A{X), whose entry Ay^i ; u, w' G V{X) 
is different from and equal to 1 if and only if v is ad- 
jacent to v' , i.e. if {vv') G E{X). It should be noticed 
that in ordinary differential geometry, graphs, like any 
other point set, can be considered as 0-dimensional real 
(or complex) manifolds. The adjacency matrix A(X) is 
related to the topology of X: a classical example is pro- 
vided by comparison between complete graphs Kg , which 
can be embedded planarly in the projective plane, but 
not in the euclidean one, and K7, which can be embed- 
ded in the 2-torus. In this sense it is possible to say that 
the graph geometry has a combinatorial origin. Never- 
theless, a graph X is still a topological space, even if not 
embedded in some differential manifold: it can indeed be 
endowed with the topology induced e.g. by the chemical 
distance metric (the chemical distance between v and v' 
being the number of edges of the shortest path connect- 
ing the two vertices) . Moreover an atlas of X over Z can 
be defined, simply by fixing an integer labeling of the 
vertices of the graph. In that case it would be possible to 
consider X as a sort of one dimensional integer manifold, 
with a non trivial topology related to its connectivity 
properties and with every vertex unambiguously identi- 
fied by one integer coordinate. 

Let us now introduce the commutative algebra A{X) 
of complex valued functions over X, together with the 
ordinary function product. A basis is provided by the set 
of vertex characteristic functions {a^jw S V{X)}, such 
that a^(u') = Sy^yi. If f{v) = fy, then one can write 
/ = J2vev{X) f-v ■ 'Orthogonality' and 'completeness' 
relations hold: Uy ■ ay' = Sy^yi ay , J2v 
where 1 is the multiplicative identity. 

Another set of generators {a(^yy'j\{vv') G E{X)} may 
be also defined, in one to one correspondence with unori- 
ented edges {vv'). A generic ' 1-form' w, belonging to the 
set n{X), will be written as t^; = w^^/^^a^j,^/^ . 

{vv')eE{X) 



€V{X) ' 
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Identification of w as a 1-form is based on the existence 
of the following structure. First three different distribu- 
tive product laws can be introduced: 



/ ^{vv') — fv ^(vv') 7 ^{vv')f fv' (^{vv') 



(1) 



together with the operator d, connecting functions in 
A{X) to their differentials, 

a{vv') = I v ^v' ' '^^ ' " ^ ■ ^ ' 

It follows from (jSJ that the generators {a(yyi) } are not lin- 
early independent, and, therefore, that different explicit 
expansions of fundamental differentials {da^ G ^{X)} 

can exist. It is evident that df = f^' cl{vv')i 

{vv')eE{X) 

whence the following antisymmetrized expansion derives: 



df — ^ ifv' - fv)a{vv') 

{vv')eE(X) 



(3) 



The two expressions differ only by a 1-form proportional 
to d{l), null by definition. One should notice, how- 
ever, that even if d{l) is the identity for the sum of 
1-forms, it is not a trivial multiplicative absorbing cl- 
ement, because the product of two null 1-forms is in 
general a 1-form different from 0. Let consider for in- 
stance = a{vv') = 1 = 0: nevertheless 

{vv')eE{X) 

one has Lu^f = ^ = iu^f — f uj-^, that is: 



wi/ = df , 



(4) 



which is not a null form unless / is constant over X. 
Thus, the infinitely many possible expansions of the dif- 
ferential are not strictly equivalent one to the other and 
we have to fix a particular way of writing it. The choice 
© is the most natural because of its manifest similarity 
with a finite differences approximation of the usual differ- 
ential. It should however be recalled that Q is an exact 
relation: the ambient space in fact is discrete and there 
is no need to take continuous limits (in other words, we 
are dealing with a generalization of the /i-calculus, see 
e.g. Quite interestingly, the definition Q yelds an 

unconventional form of the Leibnitz rule: 



difg)^ fdg + gdf + dfdg . 



(5) 



It can be easily checked that every expansion of the dif- 
ferential df produces some additional nonlinear term in 
the Leibnitz rule: this is a characteristic feature of the 
differential defined by (0), and it is not present in other 
discrete calculi known in the literature (like, for instance, 

A notion of discrete derivative can be introduced 
as well. A tangent vector field V over X, be- 
longing to the set xi^)^ is a linear combination 



V = of the generators defined by 

{vv')eE{X) 

means of the following 'duality' relation (-I-) : ^{X) x 
xiX) A{X) , (a(™/)|9(u,u,/)} = Sy^i^Sy'^w'ay . The 
action of a vector field on a scalar function is then: 

V(/)-(d/|V)= E ifv - fv)V^,,,^a, , (6) 

{vv')eE{X) 

so that i9(„t,') can be interpreted as acting as a derivative 
along the directed edge {vv'): 9(„„/)(/) = {fy' - fy)a^ . 
In a sense the derivative along a vector field V can be 
considered as a derivative along a superposition of differ- 
ent edges {vv') e E{X) with weights given by V(i,„'). 

Vector fields can also be interpreted as generators of 
transformations of A{X), turning the set x{X) into an 
algebra. Let m be a map m : A{X) —f A{X). Indeed 
setting 



m = 1 + V„ 



(7) 



one defines first Vm, referred to as the infinitesimal gen- 
erator of m, which can be assumed, with no loss of 
generality, to belong to the algebra of discrete deriva- 
tives along vector fields. The name is obviously sug- 
gested by Lie algebras theory, because of the analogy 
between Q and the first-order perturbative expansion 
of an exponential. |(7J) is nevertheless once more an exact 
relation. Consider now the endomorphisms associated 
to every single permutation tt of the vertices in V{X): 
yv G V{X) , m7r(a„) = a^(„); their group is a repre- 
sentation of the symmetric group Sp acting over the al- 
gebra A{X), where p — |Ar| is the order of the graph. 
Every has a vector field as infinitesimal generator: 
— J2v£V{x) ^tt(v),v ■ The vector fields are then 
the infinitesimal generators of the action of S'p over A{X). 
Cayley's theorem, on the other hand, states that every 
finite discrete group is isomorphic to a suitable symmet- 
ric group; to show that, one might in general resort to 
the right adjoint action. Set: 

ygei& ,letng:(&^e,ng{g')=g'g . (8) 

The {TZg\g G ©} form a representation of <& over S^0^; 
the infinitesimal generator of g being then represented by 
the vector field: 



^9 = ^■^s(i').'' 

vevix) 



(9) 



Let now focus the attention on the union of the alge- 
bra A(X) with the space of symmetric 1-forms: T*X = 
A{X) U risiX); and with the space of symmetric vec- 
tor fields: T^X — A{X) U XsiX); a 1-form or a vector 
field are symmetric if their components respect the rela- 



tions A 



{vv') 



= A 



{v'v) 



Both T,X and TtX can then be 



endowed with a symplectic structure by introducing suit- 
able Poisson brakets; for instance, in the case of T*X one 
defines the bilinear operator {■,■}: T*X x T*X A{X) 
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such that: 

{av,a,^'} = 

{a'{{vv')),aw} — {Sv'w - Swv)av (10) 
{a((j,^')),a((^^'))} = 

The symmetrization operator is a(^„^/))) — ^(a^^„/^ + 
a(„'„)). A similar bracket can be introduced in TsX 
by substituting everywhere a^„^/^ with obtaining 
then: 

{XJ} = Xif) = {clf\X) (11) 

The operators defined in H1Q|I (and (|1HI ) are actually Lie 
products, because they are antisymmetric and respect 
the Leibniz and the Jacobi relations: 

{f,iu}^-{ujj} 
{fg,u;}^{f,u;g} + {gju;} (12) 
{u;,{vJ}} + {7^,{f,u;}}=.0 

where f,g G AiX) , X e XsiX) , uj,r] e n,{X) (sym- 
metricity is a sufficient condition for proving the Jacobi 
property). It can then be stated that generators of func- 
tions and of 1-forms (vector fields) over the graph X arc 
canonically mutually conjugate; in a sense, one may claim 
that a canonical conjugation relation holds between de- 
grees of freedom related to vertices and to edges of the 
graph. 

Consider now a material particle constrained to move 
over the graph X. At each time the particle will lie in 
one vertex v G V{X) and at the following step it will 
reach one of the sites v' connected with v, such that 
{vv') G E{X) as specified by the adjacency matrix A{X). 
One can define a 'state function', say ip, belonging to 
the function space A{X), whose meaning is different ac- 
cording to the dynamical description of the system one 
aims to achieve: if the particle position can be known 
precisely at each time step, its configuration (and hence 
its state) will be simply described by a vertex charac- 
teristic function ay, which is non zero only at the site 
where the particle is located; if the particle is a classi- 
cal random walker (or is a quantum particle) the state 
function (or its square modulus) should provide the cor- 
responding localization probability density over X. On 
the other hand, the most general way to introduce a po- 
sition operator is by defining it through a parameteriza- 
tion Q : V{X) Z, , V Qy and an associated function 
Q = X]«ey(x) Iv'^v The action of the position operator 
on the state function will then simply be required to be: 

QW = Q-i^. (13) 

In particular the vertex characteristic functions will be 
eigenstates of the position operator: Q{ay) = qyUy. If 
the vertices of graph X correspond to positions eigen- 
states, its edges will have to be connected to directions 
of motion, that is to velocities or momenta: this is the 
most natural physical interpretation of the canonical con- 
jugation between functions and 1-form generators. If the 



particle hops from its initial position to a neighbouring 
site, it will cover always a single unit of chemical distance 
and every bond will be associated to the same elementary 
momentum eigenvalue, that can be set equal to 1 without 
loss of generality. An obvious choice for the momentum 
operator V will then be the symmetric 1-form uj-^ , which 
is such that cjjj^a^^,^/^ = a(«i,'): 

V{tp) = Lj^ij = dip . (14) 

The Poisson parenthesis of position and momentum oper- 
ator is then different from zero, and gives rise to canonical 
conjugation: 

{V, Q}^AQ^J2 - '?.)a. • (15) 

{vv') 

The physical meaning of the momentum operator can 
be made more trasparent if the graph X is endowed with 
a richer structure. Let X be the Cayley graph of some 
group generated by a set 9) of elements: for example, 
if <S is defined by a presentation, the set of the group 
generators provides a natural choice for io. A general 
function in A{X) will be written ip — tpg Ug and its 

gets 

differential as dip = ^^("^sft. ~ ^g) o.{g gh)- Let now 

X be a graph whose vertices are in one-to-one correspon- 
dence with the complex numbers {x{h)}, where x denote 
now the characters of the group G3 (corresponding to its 
irreducible representations) and h € ^. This is consis- 
tent because the set of generators represents the set 
of possible 'directions of motion' from each point of X: 
the vertices of the dual graph X, being labeled by the 
characters evaluated along the elements of Sj, are there- 
fore associated to possible moves over X. One should of 
course choose different generators in different conjugacy 
classes, in order to have inequivalent x(/i). On the con- 
trary, there are no constraints on E{X) (whose element 
{gg') , .9, g' G G depends on the single group element 
g" — g' g^^): nevertheless, if G is abelian, the set of its 
characters will be again a group and appropriate choices 
of the edge set will make X one of its Cayley graphs. 
The group structure allows harmonic analysis over the 
algebra A{X); a 'Fourier transform' ip — ''Px^xW of 

X{h) 

Ip can indeed be properly defined (both in the abelian 
and in the non abelian cases): 

^x^Yl ^aX* (g) , ^g"^ -^Yl X{g) , (16) 
ge& X 

where G and H are the orders of © and Sj, respectively (a 
different normalization constant will appear in the infi- 
nite group case) . Like in the case of the ordinary Fourier 
transform, Plancherel and Parseval theorems hold; given 
a function t G A{X) such that tg = ipgt, the components 
of its Fourier transform will be — x{t) ijJx'i ^^'^ for any 
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norm defined in A{X) it will be possible to introduce one 
in A{X) such that = 

Momentum can now be introduced as a sort of not in- 
jective parameterization over X, namely a linear operator 
defined in such a way that: 

V{b^[h)) ^ Px{h)b^ih) , (17) 

with x(/i) = 1 +Px{h) ■ Similarity with the customary re- 
lation e**^^ ^ l+ikx+0{k'^) should be noticed once more. 
Let now consider the case of a quantum particle. Expec- 
tation values of the square of position and of momentum 
are interesting because they give information about their 
variances: explicit calculation shows that uncertainty in- 
equalities hold in full analogy with Heisenberg Princi- 
ple. In a state ?/; — J2g£@ '^g^g^ l^ie second probability 

distribution momenta are given by {Q?)^ = \'ipg\^qg 

36® 

and by ("P^)^ = ^ j^^j^ (^(/j) _ i). One has there- 
fore: {Q^h{f^h =Y.ge0\i^g\' <llY.x,He^\^x? \x{h)- 

Once again, by Plancherel theorem |'P('0)l — 

^geejiei) |(^V'){5ft g)f = CT'(^^)'/'- ^^'^ Matter relation 
highlights the connection between the combinatorial and 
the group theoretical definition of momentum. Indeed 
{dip) (^g gii^ = fgh — fg IS thc (?-th compoucnt of the func- 
tion Afiip = d(^g g/i)'07 whose Fourier components are 

(Ahip) — Px(h)'^x'- case of an infinite group 

with infinitely many inequivalent irreducible representa- 
tions (like most Lie Groups) a one-to-one mapping is thus 
established between E{X) and V{X) and the two mo- 
mentum operators become simply proportional to each 
other. Using now general, well established inequalities 
0, one has: 

{Q^){V^) > ^ ■ (18) 

Even if H18() is not the strictest possible inequality which 
can be obtained, nevertheless it appears quite similar to 
the standard physical Heisenberg uncertainty inequality. 
This is suggestive, because we are not dealing with actual 
positions and momenta, nor with their corresponding op- 
erators as usually introduced in quantum mechanics, but 
only with their abstract analogue over X, nothing more 
than rough notions of position and direction. Neverthe- 
less this is sufficient to obtain the inequality (|18() : a hint 
of the deep ultimate origin of the fundamental Heisenberg 
uncertainty in quantum mechanics. 

The introduction of a precise notion of time evolution 
is now required for the consistent construction of a dy- 
namical system. Ambient space will be the given graph 
X, but also time will be discretized, in order to supply a 
fully covariant description. Furthermore, a discrete time 



is naturally endowed with partial ordering, provided by 
its integer labelling; in other words, there exists an in- 
trinsic set-theoretical arrow of time, associated with the 
causal ordering of events in special relativity, but also 
with the algebraic time-flows that one can construct in 
generally covariant theories (such as the thermal time of 
ref. i). 

The state function of the system will now be defined 
over a larger graph, called the time expansion t{X), ob- 
tained joining different copies of the original graph X, 
in accordance with the scheme induced by the dynamics; 
each vertex of the time expansion will be denoted by a 
double label, {v,ti), in which the first entry, v G V{X), 
will once more give the position while the other, ti, will 
have the meaning of temporal coordinate. At time ti, the 
state function will then be written as: 

Let focus, for example, on random walks and quantum 
evolutions and let assume a discrete spacetime model, 
where infinitely many edges of the type ((u, ti){v' , ii+i)), 
i G N, are introduced for each edge {vv') in E{X). The 
graph endomorphisms generated by vector fields of the 
form: 

E = X! -£^((«,t,+i)(t>',t,))5((«,t,+i)('«',t,)) > (20) 

{vv')£E(X),i 

will generically describe the system evolution, but other 
constraints must be respected if the resulting dynamics 
has to be a physically allowed one. In particular, the 
evolved state function ipi+i must be a linear combination 
only of the {0(^,^4^^^)}, because evolution is toward the fu- 
ture (in a sense, the time sequence is completely identified 
with the causal sequence itself). Thus the components 
along the past basis {a(^, j.)} have to vanish. This hap- 
pens if and only if the evolution generator E satisfies the 
progressivity condition: 

yvev{x), ^((-.t.+i)(«'.t.)> = 1 • (21) 

v'ev{x) 

Condition (|21|l guarantees that the evolved state func- 
tion is identical with its restriction to the next time- 
sheet, and there is no state self-overlap in time, 
thus implementing the right arrow of time directly 
in the evolution law: tpi+i = [1 -f E] {ipi) = 
J2v^v'ev(x) ^{{vA,+i){v',u))iPiv',u) a{y,t,+i) ■ Finally, the 
progressivity condition can be rephrased as quasi- 
stochasticity of the evolution matrix (in a quasi- 
stochastic matrix entries are in general complex and their 
sum on each row is equal to 1), whose w'-entry is equal 
to i?{(i,',t;_,.i)(„,t.)) . It is worth noticing that this con- 
straint follows only from the existence of a time flow with 
a precise direction; in the random walk case has to be 
stochastic to begin with, being the transition matrix of 
a Markov Chain. Quasi-stocasticity must subsist, how- 
ever, also in the quantum case, besides the conservation 
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of probability, associated usually to the unitarity of the 
dynamical flow. A strong connection is then evident be- 
tween quantum and diffusional dynamics, as suggested in 
the stochastic quantization literature ^Td\. 

Introduce next the temporal derivation operation 



)>■ 



The time derivative 



v£E{X) 



^(ri)ri+i) "0 can be evaluated directly, for the time- 



dependent state function ijj 



ipi of a single particle 



moving on X. If -0^ is the column vector of the compo- 
nents of ipk, it is straightforward to write: 

d{r,r,+^)i^ = «Hfc • , for Efc = 1 + aHfe . (22) 

aHfc can be interpreted as the infinitesimal generator of 
the evolution matrix E^, in analogy with what was done 
in 10. The same equation describes both the evo- 
lution of probability distributions and of single particle 
wave functions and in both cases 1 -I- aH has to be quasi- 
stochastic. For a random walk, it must be also positive, 
while for quantum evolution unitarity should be ensured. 
In this simple model quantum evolution can thus be ob- 
tained as an analytical continuation of a Markov process 
over the complex plane, that is as a generalized random 
walk with complex transition probabilities; the meaning 
of such odd entities remains far from obvious, but some 
enlightening 'physical' remarks can be found in Feyn- 
man's discussion (IL] : moreover complex valued Wiener 
measures are already known and have been rigorously 
characterized 

If one assumes now by analogy that H has the dimen- 
sions of an energy, like an actual hamiltonian operator, 
the constant a must have the form a = tq/H, where a typ- 
ical time scale (related to the duration of the elementary 
time step) has been factored out and H has the dimen- 
sions of an action as needed for dimensional consistency. 
Performing the time variable change t ^-^ —it, analogous 
to an inverse Wick rotation, (|22|l can be recast in the 
form: 



life • V'fc 



for 



1 



itn 



(23) 



Dealing with free particles, a natural choice for the ma- 
trix H might be proportional to the graph laplacian L = 
DE)-^ = JJ — A where A, D and J are respectively the ad- 
jacency and incidence matrix and the diagonal matrix of 
valencies of the ambient graph X (see for example jlij'). 
The incidence matrix provides the matrix form for the 
momentum operator 1-form oj-^ and the laplacian is then 

a sort of square of the momentum operator. H = ^L, 
is a symmetric hamiltonian; it is then possible to intro- 
duce the vector field H = Y.(vv')eE(x)^^'-"^{(v' MvA)) 
and rewrite as: 



(24) 



The constant {ih) ^ appears in front of the Poisson 
brackets only because of Wick rotation and reasons of 



dimensional consistency, and not in force of Bohr Corre- 
spondence Principle: there is no quantization to perform 
simply because the mechanics is quantum-like from the 
beginning. 

Unfortunately, 1 — ito^L, is quasi-stochastic but not 
unitary, and the evolution preserves the (normalized) 
sum of the components of tp but not its 2-norm: the state 
function can then apparently be considered a complex 
generalized probability distribution but not yet a proba- 
bility amplitude in the usual sense. Nevertheless, let us 
introduce a larger time scale t — n to, where n is a pos- 
itive integer; the hamiltonian is time-independent, then 
the evolution matrix over the full time interval t can be 
written as a Trotter-Suzuki expansion E„ (l — ^) , 
converging for n large enough to E = exp (^—j^Mt). E is 
both quasi-stochastic and unitary, then on average the 
norm of ijj is conserved and at the new time scale t the 
generalized random walker behaves exactly like an or- 
dinary quantum particle. Deviations from unitarity are 
very difhcult to observe if the original time scale to is 
very small: for instance, if to is taken of the order of the 
Planck time (10"'*'^ s), for a negligibly small renormal- 
ized time scale t ~ 10"^'' s (a 'yoctosecond') the relative 
deviation would be of the order 10~^^. At a more formal 
level the same result, hinting at the convergence of the 
dynamical matrix, comes from direct integration of 124(1 : 



■0* = exp 



^AdH 

h 



(25) 



written in term of the exponential of the adjoint action 
of the dynamical algebra. 

Quantum phenomena (oscillation, interference, local- 
ization and delocalization) are then quite easily achieved 
either as coarse-grained behaviours of a generalized diffu- 
sion or as solutions of a fundamentally classical equation 
of motion. A very interesting feature is here that tun- 
nelling between sites not directly connected by edges can 
be observed: for instance if the particle is initially local- 
ized at one of the endsites of a linear chain, after only one 
renormalized temporal step t there will be an exponen- 
tially damped but not null probability of finding it at the 
opposite lead. Of course such tunneling is the byproduct 
of a much faster sublying quasi-stochastic dynamics. The 
same happens for the related quantum probability inter- 
ference. The example problem studied in this paper is 
little more than a simple toy model and the theory is still 
far from being rubust enough as to describe real world 
physics; nevertheless sound arguments for the existence 
of a deep-lying connection between space-time discretiza- 
tion, random processes and quantum dynamics have been 
put forward. It may be reasonably hoped that quantiz- 
ing geometry may indeed help to better understand the 
very meaning of the quantization of physics. Basic pil- 
lar of such understanding is mimicking non-commutative 
geometry [Tsl l , in the context of which it was shown that 
a mathematical framework can be developed where the 
fundamental object of (continuous) geometry is no longer 
a manifold but an algebra. Since commutative C-* alge- 
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bras biuniquely correspond to locally compact topologi- 
cal spaces, all relevant information about the topological 
structure of such spaces is encoded in their algebra of 
functions. This underlying notion is what has allowed us 



to construct a 'differential calculus' over ambient spaces 
given by graphs, resorting to the differential structure 
that can be straightforwardly defined over a commuta- 
tive algebra. 
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